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Abstract 

We consider the low-energy particle-particle scattering properties in a pe¬ 
riodic simple cubic crystal. In particular, we investigate the relation be¬ 
tween the two-body scattering length and the energy shift experienced by 
the lowest-lying unbound state when this is placed in a periodic hnite box. 
We introduce a continuum model for s-wave contact interactions that respects 
the symmetry of the Brillouin zone in its regularisation and renormalisation 
procedures, and corresponds to the naive continuum limit of the Hubbard 
model. The energy shifts are found to be identical to those obtained in the 
usual spherically symmetric renormalisation scheme upon resolving an im¬ 
portant subtlety regarding the cutoff procedure. We then particularize to the 
Hubbard model, and hnd that for large hnite lattices the results are identical 
to those obtained in the continuum limit. The results reported here are valid 
in the weak, intermediate and unitary limits, and can be used for the ex¬ 
traction of scattering information ,via exact diagonalisation or Monte Carlo 
methods, of two-body systems in realistic periodic lattices. 

Keywords: Scattering theory. Effective held theory. Lattice fermions. 
Finite-size ehects 


1. Introduction 

The notion of low-energy universality in particle-particle collisions is a 
powerful concept that has been around, in one form or another, for a very long 
time. The general idea consists of replacing realistic, complicated interactions 
with much simpler ones and renormalising their bare coupling constants in 
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favour of exact (either experimental or theoretical) scattering properties. The 
model interaction, irrelevant at the two-body level, can then be used to 
investigate the effect of interactions with higher particle numbers without 
the extra complications of the realistic interactions. For example, in 1957, 
Huang and Yang introduced their s-wave regularised pseudopotential [ij, 
which aimed at reproducing the exact scattering length of a realistic two-body 
process by means of a very simple model interaction. The idea, however, was 
put forward much earlier by Fermi in the 1930s j^, who studied neutron- 
nucleon scattering, and was able to £t the scattering length in the first Born 
approximation. 

The calculation of two-body scattering properties with rather realistic 
two-body interactions, if these are spin-independent, spherically symmet¬ 
ric, and single-channel, is not too complicated, especially with the comput¬ 
ers we have access to nowadays. It is genuine multichannel collisions that 
can become quite challenging to approach directly from the corresponding 
Lippmann-Schwinger equation. In nuclear physics, there is a fundamental in¬ 
terest in obtaining nucleon-nucleon scattering properties from hrst principles 
using lattice QCD (see Q and references therein). In condensed matter and 
cold atomic physics, particle-particle scattering theory in a crystal, where the 
incident waves are Bloch waves, and with a realistic (e.g. screened Coulomb) 
interaction, constitutes a formidable problem: all bands of the periodic po¬ 
tential are coupled by the interaction, Galilean relativity does not hold so 
that the collisional properties depend on the total center of mass momentum, 
and spherical symmetry is not present and cannot be exploited. 

In situations like those described above, it is of great appeal to be able 
to extract scattering properties without having to directly use scattering 
theory. In massive quantum held theories, Luscher showed how to extract 
partial-wave scattering amplitudes, i.e. phase shifts and scattering lengths, 
from the energy shifts due to the interactions when the system is placed in 
a hnite volume j^, thereby generalising an older result from non-relativistic 
quantum mechanics [H. However, neither Liischer’s nor Huang-Yang’s results 
for the scattering lengths hold when these are unnaturally large, i.e. near 
the unitary limit. This problem was elegan^ solved using pionless effective 
held theory by Beane et al. in ref. j^. In this way, an analysis of 
low-energy nucleon-nucleon collisions, directly from QCD, was possible soon 
after [^. In the case of particle-particle scattering in a crystal, hnite-size 
ehects in the forms given in refs, and [^, however, do not appear to be 
known. In this article we derive the relevant formulas from ehective held 
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theory on a tight-binding lattice, that is, the Hnbbard model. We will focus 
on the three-dimensional (space) case, for the one-dimensional case is well 
understood jsl and Liischer’s formula in one dimension holds. We derive 
the pertinent expressions using Bethe-Goldstone equations, and we begin 
by proving these via the introduction of a more general class of methods 
consisting of kernel subtractions. We then introduce the Hubbard model 
and analyse, in first instance, its naive continuum limit. We then return to 
the Hubbard model itself and show that, at low energy and large volume, 
the energy shifts are completely analogous to those obtained from the naive 
continuum limit, which shows once more how low-energy physics, at least in 
the two-particle case, is universal. 


2. Kernel subtractions and Bethe-Goldstone theory 

For general interactions, one way of establishing the relation between the 
scattering states of a system (i.e. the positive energy states in the infinite 
size limit) and the energy shifts in a finite box is to use Bethe-Goldstone 
(BG) theory in vacuum, that is, without the static Fermi sea background. 
In fact, even if it looks very different from it, Liischer rederived, and used the 
exact same old Bethe-Goldstone’s equation in his seminal paper [^. Here we 
will relate BG theory to a much more general method - kernel subtractions - 
for solving homogeneous integral equations 0, and show that it corresponds to 
a particular choice of kernel subtraction. The method is easy to implement. 
Let \'ip) be the solution to the integral equation 

|^) = K(i?)|^), (1) 

where K{E) is the kernel (operator), with E some parameter, which in the 
relevant case here will correspond to an energy eigenvalue. The above equa¬ 
tion ([1]) is homogeneous and therefore only has a solution for particular values 
of E. Gonsider now a vector I?/;) that satisfies the following integral equation 


= \E)+ K{E)-\E){T\ 


( 2 ) 


where \E) and (F| are arbitrary vectors. Then, for the values of the energy 
for which (F|'0) = 1, it holds that \^jJ) = 


^Here, integral equation is taken in the most general sense, and applies equally to 
discrete kernels. 
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We now particularize to the physically relevant case, and derive the BG 
equation. Let satisfy the stationary Schrodinger equation 




(3) 


We assume the Hamiltonian is of the usual form H = Hq + V. Then, if E 
is not in the spectrum of Hq, the wave function satisfies the homogeneous 
Lippmann-Schwinger equation 


IV^) = Gq{E)V\^1:), (4) 

where Gq{E) = {E — Hq)~^. The kernel, Eq. ([T]), is therefore given by 
K{E) = Gq{E)V. We choose the following vectors |F) and |r), 

1^) = |k), (5) 

(r| = {k\Go{E)V, (6) 

where |k) is an eigenstate of Hq, i.e. Lfo|k) = e(k)|k). Then, |'0) satishes 

IV^) = |k) + [Gq{E)V - |k)(k|Go(E)H] IV^), (7) 

E = e{k) + {k\V\^l^). ( 8 ) 


Eq. ([7]) together with the auxiliary condition (|8]) correspond exactly to BG 
theory in operator form. In order to relate it to the finite-size reaction matrix, 
we introduce the identity 1 = X]k|k)(k| and define Vlijj) = f|k). After 
straightforward algebra, we find 


(k'|f |k) 

E 


(kWik) + 

q#k 


(k'|r|q>(q|f|k) 
E - £(q) 


e(k) + (k|f|k). 


(9) 

( 10 ) 


3. Hubbard model 

We shall focus on one specific model, namely the Hubbard model on a 
simple cubic lattice. This is the minimal model for interacting electrons (or 
any other kind of spin-1/2 fermions) in a crystal, and consists of a nearest- 
neighbour hopping, or kinetic energy term, and an on-site interaction term. 
Its Hamiltonian in second quantised form is given by 

H =-J ^ c\^Cj„ + G QJN. (11) 

(*P>,o'=t4 * 


4 



Above, J > 0 is the hopping rate, (i, j) denotes nearest neighbours, with 
i,j eZ^, Ci„ {cl) are fermionic creation and annihilation operators, U is the 
on-site interaction strength, riia = clcia is the local number operator and 
N = is the total number operator (we have set the lattice spacing 

d = 1). In Hamitonian flTT]) . the last term on the right-hand side amounts to 
an energy shift so that the single-particle ground state energy vanishes. In 
particular, the kinetic energy dispersion, e(k) of the Hubbard model flTT]) is 
given by 

e(k) = —2J [cos(/c„) — 1], (12) 

a=x,y,z 

where the quasi-momenta ka G [—tt, vr), and the Brillouin zone is given by 
BZ = [—TT, vr)^. If the system is placed on a hnite cube of side length L = Lg 
(recall the lattice spacing d = 1) with periodic boundary conditions, the 
quasi-momenta are restricted to take values ka = 27rna/L (mod 27r), with 
Ha e {0,1,... ,L - 1}. 


4. Naive continuum limit 


Before studying the Hubbard model in detail, it is very instructive to 
consider its naive continuum limit. The bare continuum Hamiltonian is given 

by 



drVV’i(r) ■ V'0a(r) +9 


/rG[0,L)3 


drq(r)q(r)Vii(r)Vif(r). 


(13) 

Above, m is the effective mass, g is the bare interaction strength, and 'ipa 
{'ipl) is the fermionic annihilation (creation) operator in the continuum. 

The above Hamiltonian, Eq. flT^ . is obtained as follows. Firstly, we 
restore the lattice spacing d (which was set to d = 1 in the previous section), 
in such a way that the action of the non-interacting Hamiltonian on a single¬ 
particle wave function, in the hrst quantisation, is given by 


{Ho'ijj){r) =-J ^ [ip{r + ead) + 'ijj{r-ead)-2'ijj{r)], (14) 

a=x^y,z 


where r = {x,y,z) = {nx,ny,nz)d, with Ua G Z, and where we have dehned 
the unit vectors Cq such that Cq ■ r = a {a = x,y,z). The continuum limit 
of Eq. flTT)) is attained by setting J = /2md?. The continuum limit of the 

bare interaction strength is attained by identifying U = g/d^ as d ^ 0. 
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The effective interaction in Eq. fll3p is a zero-range ”s-wave” interaction. 
It is renornialisable, and accounts correctly for the s-wave scattering length 
with a vanishing effective range. This interaction can be renormalised in a 
number of ways, all equivalent to each other. For instance, minimal subtrac¬ 
tion in the hard cutoff regularisation scheme [sj, dimensional regularisation, 
Huan g’s pseudopotential [l|, or the more recent method of Tan’s distribu¬ 


tions [id, llll , all give the same results. In all the aforementioned methods, 
spherical symmetry of the contact interaction is exploited. The question 
arises of whether the cubic symmetry of the lattice Brillouin zone in the 
continuum limit changes the status of affairs. The answer to this question, 
which seems quite trivial, is not straightforward at all. On the one hand, 
the scattering properties at low energies remain unchanged. On the other 
hand, Liischer’s formula heavily relies on spherical symmetry, a fact which 
is particularly evident in the pionless EFT considered by Beane et al. in ref. 
[^, and we will show that Luscher’s formula does change significantly when 
the lattice symmetry is respected in the regularisation process, unless special 
care is taken in the renormalisation process. 


4-1- Scattering theory 

In this subsection, we shall obtain the reaction matrix (i?-matrix), by 
regularising the contact interaction in such a way that the symmetry of the 
Brillouin zone is manifestly preserved at all steps. The regularised contact in¬ 
teraction V (k, k') in the momentum representation for a spin-singlet fermion 
pair with ultraviolet (UV) cutoff A is therefore given by 

(k'|E|k) ^ E(k, k') = ge{A, k')0(A, k), (15) 


where we have defined 


9(A.k)= n (16) 

a=x,y,z 


Above, 9{q) is the Heaviside step function. The above interaction, Eq. (IT^ . 
has to be contrasted with the widely used spherically symmetric regularisa¬ 
tion 0,0 [3, which reads g6{A? — k^). 

After separation of center of mass and relative coordinates, and at van¬ 
ishing center of mass momentum, the bare i?-matrix i?*with the interaction 
in Eq. flT^ is straightforward to calculate, and reads 


(k'|.R,(E)|k) 


0(A,k')0(A,k) 

l/g + h{E)/{2nr 


(17) 
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where 


Io{E) 


JqG[-A,A)3 


dq 


E — fi?q^ /m' 


(18) 


with V denoting Canchy’s principal valne. The UV strnctnre of Io{E), Eq. 
dUD, can be extracted from its valne at zero energy. Nnmerically, we have 
obtained /o(0) = —15.3482484734169.. .A + C>(l/A). The zero-energy on- 
shell i?-matrix can be renormalised in favour of the scattering length a by 
means of minimal snbtraction, i.e. by setting 


1 _ 1 /o(0) 

9 9r (2vr)3’ 


(19) 


in which case the renormalised conpling constant takes the valne gji = 
Anh'^a/m. The i?-matrix at arbitrary energy is renormalisable if Iq{E) — Io{0) 
is finite for all energies E, that is, if the linearly divergent part of the integrals 
Io{E) is energy-independent. To show that this is indeed the case, write the 
integral Io{E) as 


Io{E) 


V 



dq 


1 

E — h?q‘^/m 


'qG[-A,A)3 


6{q‘^ — A^) 
E — h?q‘^/m 


( 20 ) 


The second integral on the right hand side of Eq. (1201) can be calcnlated as 
a series expansion that is convergent for any finite energy. We expand the 
integrand as 


e{q^ - A2) 0(g2 _ A2) 






mE 


E — h?q‘^/m 

The resnlting integral in Eq. (|2ni) is therefore given by 


( 21 ) 


I„(E) = /„(0) 


m 

¥ 


E 


dqe{q^ - A^ 


k 


2n 


= 1 JqG[-A,A)3 


2(n+l) ’ 


9 


( 22 ) 


where /c^ = mE/h^. The integral of the n-th term in the series is bonnded 
as 


0 < 


dqe{q^ - A^ 


< 


47r 


Jqe[-A,Ad (2n-l)A2’^-i' 

Therefore, the following bonnd holds 


Jcie[-A,A)^ ^^E - K^q^/m 


0(g2 - A2) , . 


n=l 


(k/A)‘ 


(23) 


(24) 
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In the limit A —)■ cxo, the integral is bonnded by zero from above and below, 
and it therefore converges to zero. Conseqnently, Io{E) = /o(0) and the 
renormalised i?-matrix in this regnlarisation scheme, for finite energies, is 
constant and identical to that obtained by means of spherically symmetric 
regnlarisation. 

There is a subtlety that needs to be considered in order to obtain mean¬ 
ingful results. The integer cutoff in lattice sums with Lg lattice sites (which 
we consider odd for convenience) is A = {Lg — l)/2, which implies that the 
lattice site number-dependent real cutoff A{Lg) is given by 


A{Lg) 


71 

d 



(25) 


The interaction must be renormalised by using the L ^ oo limit of the cutoff, 
A(cx 3), since the limit L —)■ oo is taken at the start of the calculation when 
doing scattering theory, that is, the renormalisation prescription reads 

1 1 1 Tib m 1 m/h^ , 

where b = —15.3482484734163 .... Obviously, in the limit L —)■ oo, we recover 
the result from minimal subtraction. 



4-2. Liischer’s formula 

We place the system in a cubic box of side length L with periodic bound¬ 
ary conditions, at vanishing pair momentum K = 0. We will consider here 
the lowest unbound state. From BG theory, Eq. (|9]) and Eq. (HB, and 
setting E = e(0) -|- AE = AE, we obtain 


AE 


A + A 

L3 L3 






AE — e(q) ’ 


(27) 


where the sum is assumed to be regularised and is renormalised to any given 
order by using the renormalisation prescription for the bare coupling constant 
(?, Eq. (j26]). To next-to-leading order in g/L^, the energy shift reads 

" h - ip) 


\-in • 


where we have defined 


Sn(L) = 


1 


(29) 






In the above sums, Si is regularised with a cutoff A{Ls) (given by Eq. (12^ . 
and will be seen to be renormalised by subtracting /o(0)), calculated using 
the cutoff A(cx)), while for u > 1 the sums Sn are regular and convergent. If 
we write 


^i(^) 

L3 


m 

(27r)^h2 


/ ^^4+ 



(30) 


we obtain, to lowest order 


S'i(L) m 
L3 ~ h2(27r)2 



where a ~ —1.23951. Upon renormalisation, we obtain 


(31) 


AE 


47rh^ a 
m 


1 + 2.837297y 

1j ■ 


(32) 


The above Equation fl32|) is in perfect agreement with Liischer’s formula [^, 
and emphasizes the importance of the subtlety regarding the renormalisation 
of the effective interaction in the naive continuum limit, as discussed in the 
previous section. 

There are several advantages of working with Eq. fl3^ with respect to the 
way one obtains Liischer’s formula in the continuum limit of a cubic lattice 
[^. Firstly, there is no need to project the interaction and the i7-matrix onto 
their s-wave components, and since Eq. fl5^ only depends on the scattering 
length, this can be extracted with ease. Secondly, ’’brute force” summations 
(i.e. direct computation of the sums, equivalent to exact diagonalisation) 
provide smooth convergence for all the sums Sn{L) in Eq. (|29|) . and equally 
for the energy-dependent sums in Eq. (l27|) . This is to be contrasted with 
the direct computation of the sum analogous to Si (L) in the s-wave channel 
(3, [^. In that case, the sums are restricted to q < A, while in the present 
case these are restricted to \qa\ < A (a = x, y, z). In Fig. [T]we plot the value 
of the sums 


n^O 


0(A,n) 


-&A, 


Luescher 

'^1 


w=i: 


n^O 


e{x^ - 

n 

P) 

n ' 

2 


(33) 

(34) 


as functions of the cutoff A, where we have subtracted A = A(L) = A(cxo) — 
1/2, while the second sum is renormalised with the spherically-symmetric 
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Figure 1: Left: Blue dots are the numerically computed sums gLuescher g^g ^ function of the 
integer cutoff A (see Eq. and the dashed red line is its asymptotic value (A ^ oo) 

quoted in refs. [J,l6j. Right: Blue dots are the numerically computed sums si as a function 
of A (see Eq. (|33|) 1. and the dashed red line is the extrapolated value for A —oo. 


prescription of ref. j^. Notice that the standard result for Liischer’s sum 
is obtained here as —1.23951 — 6/2 = —8.91363.... As is clearly observed 
in Fig. [H direct computation of the sums using projections of the interac¬ 
tions onto the s-wave channel gives rise to potentially high numerical errors. 
These would be inevitable in any computation relying on exact diagonalisa- 
tion techniques, and would result in poor estimates of the scattering length. 


4-3. Unitary limit 

In the unitary limit, i.e. for a —?■ ±oo, Eq. fl3^ is obviously incorrect. 
In the spherically symmetric case, Beane et al. (6| used pionless effective 
held theory to derive the hnite-size correction to the ground state energy 
when the scattering length is unnaturally large. This corresponds to the 
limit a/Li^ —)■ ±cxo. Here, we derive the corresponding expressions for the 
continuum limit of the lattice model in two different cases: (i) Beane et aUs 
case {a/L^ ±cxo) and (ii) the more natural case of a = 0{L) in the limit 
L ^ oo. 

For unnaturally large scattering lengths, Eq. flTr|) gives 


AE 


y —— 

^AB-qq) 


(35) 


Clearly, we have AE = 0{1/L‘^) in this case, and therefore AE is of the 
same order as the dispersion e(q) at low momenta. The energy shift in this 
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case is given by 


[0.09590 ... +C>(l/aL)]. 


(36) 




Again, the leading order term is identical to the one obtained in the spheri¬ 
cally symmetric case j^. 

In the case of large scattering lengths that scale as the length of the 
system, a = ckL, with a = C>(1), Eq. (1271) gives 


^E 


Airffa 1 


(37) 


Once more, AE = and by setting AE = (27r)^/3/L^, Eq. 

becomes 

13 = --^-. 

vrl- (a/vr) En^o(/^ “ 

The weak-coupling expansion in this case has the form 


m 

(38) 


2.837297a]. 

TT 


(39) 


Comparing Eqs. (15^ and (159]) . it is easy to observe that it is not always 
possible to distinguish the cases (i) and (ii) referred to above, corresponding 
to unnaturally large and ’’naturally” large scattering lengths, respectively, 
based on only one box size in the lowest energy shift. It can be tested, how¬ 
ever, by choosing two different lattice sizes without changing the interaction 
parameters. 


5. Lattice case 

In the previous section we introduced the analog of Liischer’s formula in 
the weak-coupling and unitary limits for the naive continuum limit 
of the Hubbard model. In this section we study these expressions in the 
Hubbard model itself. 

5.1. Scattering theory at low energy 

The Hubbard model, Eq. flTTl) . is a one-parameter theory, since the 
physics of the model only depends on the ratio U/J. At zero energy, its 
i?-matrix gives the scattering length. It does, however, have non-zero shape 
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parameters and, of course, higher (even) partial waves, too. On the other 
hand, since it only depends on 17/ J all shape parameters depend exclusively 
on the scattering length [^, and so do the energies. In this way, scattering 
lengths can be extracted from energy shifts in finite lattices. 

The on-shell zero-energy 72-matrix is trivially shown to take the value 


R 


1 

i/u -wi^y 


(40) 


where 


W(E) = [ dci ^ (41) 

(27r)'^ Jbz E - eo(q) 

Above, the energy dispersion becomes that of the Hubbard model, Eq. (IT^ . 
for the relative motion of two particles at zero total momentum, i.e. 


eo(q) = -4J [cos(g„) - 1]. 


(42) 


a=x,y,z 
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The zero-energy integral is known exactly and was computed by Watson 
It has the numerical value JIT(O) = —0.12636655... = uo- Therefore, the 
scattering length is given by 


m 

h?a 



W{0) , 


(43) 


where m is the effective mass of the particles, related to the tunnelling rate 
by J = ty/2m (d = 1). 


5.2. Energy shifts 

The energy shift of the lowest unbound state is readily obtained from BG 
theory, Eq. ([9]) and Eq. oa. and reads in this case 


AE 


U/L^ 

1 _ T. V 1 

L3 AE-eo{q) 


(44) 


In terms of the scattering length, Eq. fl43|) . Eq. (H4|) becomes 


AE 


Airhfa 

mL^ 


1 -|- 4:7ia{2uo) 


Anh/Cb V 1 

mL^ ^ AE — eo(q) 

q^O 


( 46 ) 


^This is analogous to what happens with a hard-sphere interaction 


fl. 
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In the weak-coupling limit, we therefore obtain 


AE 


mL^ 


1 — 47ra(2ci;o) + 


47rh^a 


E 


mL^ ^ AE 

q,j^0 


«o(q) 


(46) 


The sum on the right hand side of Eq. fl46|) can be estimated as follows. 
Consider a cubic region hi = [—ri,ri]^ — {(0,0,0)}, with > 0 being much 
smaller than the Brillouin zone’s extent {rj ir). The sum can be split into 
the two regions BZ — hi and hi. The sum over BZ — hi can be approximated 
by an integral, and therefore 

1 1 I r .1 1 1 

L3 


E 

q^O 


eo(q) (27r)= 


dq 


'Bz-n 


1 1 


E 


(47) 


In the region hi, the full lattice dispersion can be approximated by 

eo(q) 2Jq^ = 2J|n|^ 


setting 7] = {2nN/L), we obtain in the limit N ^ oo with N/L ^ 0 

^ I 

27r 


E 


eo(q) 


lV 1 ^e{N,ii) f lY 1 . - . 


n 


(48) 


(49) 


qSO ^^ / n^O 

where b = 15.3482484734169... and Si = —1.23945..., as calculated in the 
previous section. The integral in the region BZ — hi is easily calculated as 

2 


/BZ-o^^eo(q) 


/Bz^^eo(q) 


2J 


'n 


dq-^ = 

qz 


'BZ 


eo(q) \27rJ 


/£_Y b{N + 1/2) 


Using the above results in Eq. fl46]) . the energy shift becomes 

47rh^ a 


AE 


m 


L3 L 


1 + 2.837297- 


2J 

(60) 


(51) 


The above expression coincides with the corresponding formula in the naive 
continuum limit, Eq. fl32|) . as expected. However, Eq. (1^ has been obtained 
without assuming the continuum limit at all energies, and it proves that the 
naive continuum limit - with its associated subtleties - provides an adequate 
description of low-energy scattering in the weak-coupling limit. 

The corresponding expressions in the unitary limits (both the natural and 
unnatural cases) are obtained in a fashion completely analogous to the one 
leading to Eq. (|5T|) and are identical to the expressions obtained using the 
naive continuum limit in the previous section. 
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6. Conclusions and outlook 


Besides purely theoretical implications, the results reported in this work 
have interesting applications. For instance, two-particle scattering in a real¬ 
istic three-dimensional lattice is of great relevance for condensed matter and 
cold atomic systems in optical lattices, as it is necessary in order to renor¬ 
malise the coupling constant {U / J) in the Hubbard model. Unfortunately, 
the calculation of low-energy scattering properties of this system via direct 
use of collision theory is generally very involved, the reason being the genuine 
multi-channel structure of these systems. One way to estimate the scatter¬ 
ing properties is to solve for the ground state energy of the realistic model 
consisting of only a few lattice wells (the number of wells being L^). With 
the formulas derived in this article, one can ht the scattering length, and 
therefore the interaction coupling constant, to match the energy shift in the 
realistic calculation in a small box. Generalisations of the current work to 
study, e.g. effective interparticle interactions in graphene [l^ and graphene¬ 
like optical lattices [l^, are possible by using the methods presented here. In 
intrinsic graphene, as opposed to the cold atomic case and extrinsic graphene 
(lij i , however, the inclusion of the Coulomb interaction is necessary and the 
problem is more complicated, but can be done by extending the methods 
for continuous theories [l^. In particular, the rigorous construction of the 
continuum quasi-relativistic (it is not Lorentz invariant) effective held theory 
around the Dirac points is still an open problem [^. 

It will be interesting to study the so-called Bertsch parameter 1^, dehned 
as the ratio between the ground-state energy of a resonantly interacting three- 
dimensional Fermi gas in the s-wave channel to the non-interacting ground 
state energy, but for the Hubbard model itself. This can be done using 
quantum Monte Carlo simulations 1^, 12^. Given that the energy shift. 


considered here, for two fermions at the lattice resonance is the same as that 
obtained from the pure s-wave results, it is to be expected that the lattice 
version of Bertsch parameter will be identical, or similar, with its continuum, 
pure s-wave counterpart. Its calculation would actually be much less involved 
since the elimination of higher-partial waves and effective range, cleverly done 


in ref. [19|, is not necessary at all. 


^The two-body problem in the naive continuum limit has been studied in, e.g. ref. 
There, however, the interaction is purely phenomenological. 
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